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Invasion  of  hepatocytes  by  Plasmodium  sporozoites  depos¬ 
ited  by  Anopheles  mosquitoes,  and  their  subsequent  transfor¬ 
mation  into  infective  merozoites  is  an  obligatory  step  in  the  ini¬ 
tiation  of  malaria.  Interactions  between  the  sporozoites  and 
hepatocytes  lead  to  a  distinct,  complex  and  coordinated  cellular 
and  systemic  host  response.  Little  is  known  about  host  liver  cell 
response  to  sporozoite  invasion,  or  whether  it  is  primarily 
adaptive  for  the  parasite,  for  the  host,  or  for  both.  Our  present 
study  used  gene  expression  profiling  of  human  HepG2-A16  liver 
cells  infected  with  Plasmodium  falciparum  sporozoites  to 
understand  the  host  early  cellular  events  and  factors  influencing 
parasite  infectivity  and  sporozoite  development.  Our  results 
show  that  as  early  as  30  min  following  wild-type,  non-irradiated 
sporozoite  exposure,  the  expressions  of  at  least  742  genes  was 
selectively  altered.  These  genes  regulate  diverse  biological  func¬ 
tions,  such  as  immune  processes,  cell  adhesion  and  communica¬ 
tions,  metabolism  pathways,  cell  cycle  regulation,  and  signal 
transduction.  These  functions  reflect  cellular  events  consistent 
with  initial  host  cell  defense  responses,  as  well  as  alterations  in 
host  cells  to  sustain  sporozoites  growth  and  survival.  Irradiated 
sporozoites  gave  very  similar  gene  expression  pattern  changes, 
but  direct  comparative  analysis  between  liver  gene  expression 
profiles  caused  by  irradiated  and  non-irradiated  sporozoites 
identified  29  genes,  including  glypican-3,  that  were  specifically 
up- regulated  only  in  irradiated  sporozoites.  Elucidating  the  role 
of  this  subset  of  genes  may  help  identify  the  molecular  basis  for 
the  irradiated  sporozoites  inability  to  develop  intrahepatically, 
and  their  usefulness  as  an  immunogen  for  developing  protective 
immunity  against  pre-erythrocytic  stage  malaria. 
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Malaria  transmission  occurs  when  Plasmodium  sporozoites 
from  the  salivary  glands  of  female  Anopheles  mosquitoes  are 
inoculated  into  vertebrate  hosts  during  a  blood  meal.  Sporozo¬ 
ites  quickly  reach  the  liver  through  the  circulation  and  traverse 
through  several  hepatocyte  cells  by  membrane  disruption 
before  invading  and  settling  down  in  a  final  hepatocyte  for  their 
liver  stage  development.  Exo- erythrocytic  forms  of  the  parasite 
grow  within  hepatocytes  to  produce  several  thousand  merozo¬ 
ites,  which  exit  the  infected  hepatocytes  and  invade  erythro¬ 
cytes  to  initiate  clinical  malaria. 

Infection  by  a  pathogen  triggers  a  complex  and  distinct  set  of 
cellular  and  systemic  events,  some  of  which  may  be  orches¬ 
trated  by  the  parasite  to  support  development,  and  others  rep¬ 
resenting  a  host  defense  response.  Interactions  between  host 
and  pathogens  are  diverse  and  are  regulated  in  specific  patterns 
by  unique  molecules  and  mechanisms  involving  activation  of 
transcriptional  events  of  innate  and  adaptive  immunity  (1-4). 
In  malaria,  this  complex  interaction  between  pathogen  and 
host  is  a  critical  factor  in  determining  the  progression  and  out¬ 
come  of  the  development  in  liver  of  the  parasite.  Understanding 
global  changes  that  occur  both  at  the  host  and  parasite  tran- 
scriptome  level  will  allow  better  understanding  of  the  various 
host  and  parasite  factors  influencing  infectivity  and  parasite 
development  in  mammalian  hosts.  Little  is  known  about  host 
liver  cell  response  to  sporozoite  invasion,  although  consider¬ 
able  progress  have  been  made  in  elucidating  the  parasite  genes 
and  proteins  that  are  important  for  liver  infection  (5-13).  The 
few  other  studies  on  this  area  so  far  were  focused  on  specific 
factors,  such  as  hepatocyte  growth  factor  (14)  and  CD81  (15), 
whereas  studies  using  genome-wide  microarray  approaches 
were  predominantly  performed  using  non-human  host/sporo- 
zoite  pairings.  Because  the  biology  of  non-human  host/sporo¬ 
zoite  combinations  differs  from  those  of  human,  the  relevance 
of  the  findings  from  these  studies  remain  to  be  confirmed. 

The  development  of  an  effective  malaria  vaccine  has  been  an 
area  of  intense  research  in  recent  years.  Attenuated  whole 
sporozoites,  their  constituent  proteins,  or  the  genes  encoding 
them  are  useful  candidates  for  development  of  a  malaria  vac¬ 
cine.  These  studies  were  aimed  at  discovering  specific  potential 
targets  for  development  of  a  pre-erythrocytic  vaccine  or  drug 
for  malaria.  Limited  studies  have  been  conducted  on  changes  in 
hepatocyte  transcriptome  or  proteome  on  the  initial  encounter 
with  malaria  sporozoites.  The  liver  is  a  critical  gateway  for 
invading  sporozoites  that  can  be  targeted  for  intervention 
before  the  development  of  malaria  sporozoites  into  disease  pro- 
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ducing  blood-stage  parasites.  The  various  hepatic  factors  influ¬ 
encing  binding  and  invasion  of  sporozoites  into  hepatocytes 
have  being  studied  in  mice  using  Plasmodium  yoelii  and  Plas¬ 
modium  berghei  sporozoites,  but  the  biology  of  these  sporozo¬ 
ites  is  different  from  human  malaria  parasites,  such  as  Plasmo¬ 
dium  falciparum  or  Plasmodium  vivax . 

We  used  a  cell  culture  model  of  HepG2-A16  hepatoma  cells 
(16,  17)  infected  with  isolated  P.  falciparum  sporozoites,  and 
evaluated  global  hepatic  gene  expression  changes  by  high-den¬ 
sity  microarray  profiling.  Microarrays  can  identify  genome¬ 
wide  transcriptional  events  that  underlie  liver  responses  to 
sporozoite  invasion,  and  provide  insights  into  molecular  events 
and  hepatic  factors  related  to  sporozoite  invasion  and  develop¬ 
ment.  The  in  vitro  liver  cell  culture  system  allows  evaluation  of 
much  earlier  and  well  defined  stages  of  sporozoite  invasion  that 
are  not  possible  in  human  subject  studies.  Although  genome¬ 
wide  microarray  analyses  have  been  used  to  study  malaria 
sporozoite-hepatocyte  interactions  (18-20),  these  studies  are 
limited  to  either  transcriptome  profiling  of  the  sporozoites  or  of 
non-human  host  cells/parasite  combinations.  Our  present 
study  was  aimed  at  specifically  identifying  changes  in  human 
hepatocyte  transcripts  differentially  regulated  following  early/5. 
falciparum  sporozoites  infection.  We  therefore  used  the 
HepG2-A16  human  hepatoma  cells  that  have  been  extensively 
used  in  P.  falciparum  sporozoite  invasion  assay  studies  (21-24) 
and  to  which  P .  falciparum  sporozoites  are  known  to  attach  and 
invade  without  further  exoerythrocytic  stage  development.  Dif¬ 
ferential  effects  on  liver  gene  expression  profiles  caused  by 
wild-type  and  radiation-attenuated  sporozoites  were  also 
determined. 

EXPERIMENTAL  PROCEDURES 

Preparation  of  Malaria  Sporozoites — Sporozoites  were  iso¬ 
lated  from  the  same  batch  of  Anopheles  stephensi  mosquitoes 
infected  with  either  irradiated  or  non-irradiated  P.  falciparum 
(NF  54  strain)  using  a  modified  microcentrifugation  technique 
described  previously  (25).  Briefly,  mosquitoes  were  placed  on  a 
glass  slide,  head  and  thoraxes  were  separated,  and  the  abdo¬ 
mens  discarded.  Up  to  50  heads  and  thoraxes  were  placed  on  a 
0.5-ml  microcentrifuge  tube  previously  perforated  with  a 
25-gauge  needle  and  plugged  with  glass  wool.  This  tube  was 
placed  in  a  1.5-ml  microcentrifuge  tube,  and  50  fii  of  buffer 
(PBS,  pH  8, 1.5%  glucose,  5%  fetal  bovine  serum)  was  added  and 
the  tubes  were  centrifuged  for  2  min  at  10,000  X  g  in  an  Eppen- 
dorf  microcentrifuge.  Filtrates  were  collected,  and  centrifuga¬ 
tion  was  repeated  after  the  addition  of  another  50  p\  of  buffer. 
Pooled  filtrate  was  adjusted  to  1  ml  with  buffer  and  layered  on 
top  of  a  DEAE-cellulose  (DE-52,  Whatman)  column  to  remove 
mosquito  debris  and  microbial  contaminants  (26).  The  column 
was  prepared  using  a  25-ml  pipette  plugged  with  glass  wool  and 
packed  with  10  ml  of  a  50%  DEAE  slurry.  Sporozoites  were 
eluted  from  the  DEAE  column  by  gravity  using  20  ml  of  buffer, 
and  the  initial  14-15  ml  of  eluate  was  collected  and  centrifuged 
to  pellet  the  sporozoites.  Isolated  sporozoites  were  resus¬ 
pended  in  DMEM  and  counted  using  a  hemocytometer. 
Because  there  is  the  possibility  that  residual  contaminating  sal¬ 
ivary  gland  proteins  in  the  sporozoite  preparations  may  con¬ 
tribute  to  differential  gene  expressions,  a  mock  sporozoite 
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FIGURE  1 .  Experimental  design  to  determine  the  early  liver  responses  to 
malarial  sporozoite  infection.  Cultured  HepG2-A16  cells  were  treated  with 
wild-type  sporozoites  ( Spor },  radiation-attenuated  sporozoites  ( irradSpor ),  or 
non-sporozoite  containing  salivary  gland  extract  ( SalGId)  for  30,  90,  and  180 
min,  following  which  cells  were  washed  and  lysed  for  RNA  extraction.  Non- 
treated  cells  ( Control )  were  also  included  as  comparison. 

preparation  from  non-infected  but  the  same  batch  of  mosqui¬ 
toes  was  also  made  using  similar  scaled-down  procedures.  This 
non-sporozoite  containing  salivary  gland  extract  was  used  at 
the  same  dilution  as  the  sporozoite  preparations. 

Irradiation  of  Malaria  Sporozoites— Anopheles  stephensi 
mosquitoes  infected  with  P.  falciparum  (NF54  strain)  sporozo¬ 
ites  were  exposed  to  150  gray  dose  of  y-radiation  using  a  60Co 
source.  The  time  of  exposure  of  infected  mosquitoes  to  achieve 
the  target  radiation  doses  was  based  on  the  calibration  of  the 
irradiator  by  dosimetry  and  the  half-life  of  60Co. 

Infection  ofHepG2-A16  Cells  in  Vitro— HepG2-A16  human 
liver  cells  were  maintained  in  complete  Dulbecco’s  modified 
Eagles  medium  (DMEM)  buffered  with  bicarbonate  and  sup¬ 
plemented  with  10%  fetal  bovine  serum  at  37  °C  in  a  5%  C02 
environment.  For  infection  by  P.  falciparum  sporozoites,  1  X 
106  cells  were  cultured  in  DMEM  in  6-well  plates  at  80%  con¬ 
fluence  and  exposed  to  freshly  prepared  1  X  106  sporozoites. 
The  use  of  equivalent  numbers  of  sporozoites  to  HepG2-A16 
cells  is  to  ensure  uniformity  and  consistency  of  infectivity 
across  all  samples.  Cells  were  exposed  to  sporozoites  for  differ¬ 
ent  time  periods  (30,  90,  and  180  min),  after  which  the  sporo¬ 
zoite-containing  medium  was  aspirated  from  each  well.  Cells 
were  washed  three  times  with  complete  DMEM  medium  to 
ensure  removal  of  free  sporozoites  on  the  cell  surfaces,  and 
harvested  for  RNA  extraction.  For  controls,  uninfected  cell  cul¬ 
tures,  and  cell  cultures  exposed  to  salivary  gland  extract  pre¬ 
pared  from  non-infected  mosquitoes  were  used.  All  infection 
experiments  were  performed  in  duplicate;  however,  cells  from 
these  two  experiments  were  later  pooled  together  prior  to  RNA 
extraction  to  obtain  sufficient  hepatic  RNAs  for  microarray 
analysis.  In  each  experiment  there  were  4  sample  groups  (Fig. 
1);  wild-type  sporozoite-infected  HepG2-A16  liver  cells  {Spor), 
irradiated  sporozoite-infected  liver  cells  (. IrradSpor ),  non-in¬ 
fected  control  liver  cells  {Control),  and  liver  cells  exposed  to 
salivary  gland  extract  from  uninfected  mosquitoes  {SalGId). 
The  90-min  pooled  sample  for  IrradSpor  was  lost  on  storage, 
and  was  not  included  in  the  microarray  analysis. 

RNA  Isolation,  Amplification,  Labeling  and  Hybridization — 
Following  infection,  cell  cultures  were  washed  with  PBS,  lysed 
in  4.5  M  guanidine-hydrochloride  lysis  buffer,  and  the  total 
RNA  extracted  using  High  Pure™  RNA  isolation  kit  (Roche 
Applied  Science).  For  precaution  against  possible  RNA  deg¬ 
radation,  all  extracted  RNAs  were  treated  with  50  units  of 
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FIGURE  2.  Clustering  analysis  of  951  DEGsfrom  1 1  HepG2-A1 6  cell  culture  samples  representing  four  study  groups.  A,  hierarchical  clustering  of  samples 
using  clustered  correlation  with  average  linkage.  6,  principal  component  analysis  of  samples;  Control  (A),  non-infected;  Spor(*),  wild-type  sporozoite  infected; 
IrradSpor  (•),  irradiated  sporozote  infected;  SalGId  (■),  mosquito  salivary  gland  extract  exposed. 


ribonuclease  inhibitor  (RNasin,  Promega,  Madison,  WI)  and 
stored  at  —  80  °C.  Extracted  RNAs  were  quantitated  by 
absorbance  measurements  at  260  nm  using  the  NanoDrop 
spectrophotometer  (NanoDrop  Technologies,  Wilmington, 
DE).  RNA  fidelity  was  verified  using  the  Agilent  2100  Bio¬ 
analyzer  (Agilent  Technologies,  Santa  Clara,  CA)  and  graded 
using  RNA  integrity  number  to  assess  integrity  and  ensure 
consistency  across  all  samples.  Only  RNAs  that  conformed 
to  the  recommended  quality  control  metrics  were  used  for 
probe  preparation  and  array  hybridization.  Probe  labeling 
was  performed  using  standard  Affymetrix  approved 
RT-PCR,  in  vitro  transcription,  and  fragmentation  methods. 
Array  hybridization  was  performed  using  instrumentation 
consistent  with  the  Affymetrix  Integrated  GeneChip  System, 
and  the  Affymetrix  GeneChip  Operating  System  software 
was  used  for  probe  intensity  measurements.  Probe  labeling, 
array  hybridization,  and  intensity  measurements  were  per¬ 
formed  by  Asuragen  (Austin,  TX). 

Gene  expression  profiling  was  made  using  Affymetrix 
Human  Genome  133  Plus  2.0  GeneChip  array.  This  microar¬ 
ray  contains  over  54,000  probe  sets  representing  over  47,000 
transcripts  and  variants,  including  38,500  well  characterized 
human  genes. 

Analysis  of  Microarray  Date— Intensity  signals  from  the 
hybridized  GeneChips  were  analyzed  using  Affymetrix  MAS 
5.0  software  to  generate  CEL  (fluorescence  intensity)  files.  The 
quality  of  hybridization  signals  was  assessed  using  the  array 
outlier  percentage  as  determined  by  the  dChip  (version  1.3) 
software  (27).  In  addition,  microarray  signals  were  evaluated  for 
the  presence  of  60  liver-specific  genes  (Gene  Enrichment  Pro¬ 
filer)  as  well  as  six  short  half-life  RNA  transcripts  (SMAD7, 
IL10RA,  CXCR4,  GAT  A3 >  MYC,  RBI)  to  ensure  that  the 
microarray  signals  were  not  compromised  by  selective  loss  of 
liver-specific  transcripts.  CEL  files  were  normalized  at  the 
probe  level  by  using  the  robust  multichip  average  method  (28). 
This  method  fits  a  robust  linear  model  to  the  probe  level  data, 
analyzing  each  hybridized  chip  within  the  context  of  other 
chips  in  the  experiment.  The  average  intensity  of  each  probe 
was  expressed  as  log2.  Statistical  analysis  of  gene  expression 
data  as  well  as  gene  set  comparison  analysis  were  performed 


using  the  R  statistics-based  BRB  Array  Tools  (version  3.8) 
developed  by  Dr.  Richard  Simon  and  Amy  Peng  Lam  (29),  and 
GeneSpring  software  (Agilent  Technologies,  CA). 

RESULTS 

Effects  of  Treatment  and  Post-infection  Time  on  HepG2-A16 
Gene  Expression — Initial  filtering  and  subsetting  of  gene 
expression  data  were  performed  based  on  a  1.5-fold  change 
from  the  median  in  at  least  25%  of  the  arrays  and  showing  a  log 
intensity  variation  of  p  <  0.001.  A  total  of  1570  probe  sets, 
which  represents  1115  genes,  passed  the  filtering  criteria  and 
were  retained  for  subsequent  statistical  analysis.  Regression 
analysis  of  the  time  series  (i.e.  30,  90,  and  180  min)  expres¬ 
sion  data  among  the  four  liver  cell  culture  groups  ( Control , 
Spor  (wild-type  sporozoites),  IrradSpor  (radiation-attenu¬ 
ated  sporozoites),  and  SalGId  (salivary  gland))  indicated  that 
effects  of  sporozoite  exposure  time  on  gene  expression 
changes  within  each  group  were  minimal,  and  that  the 
effects  of  different  treatments  (i.e.  Control,  SalGId ,  Spor ,  and 
IrradSpor)  were  the  predominant  contributor  to  gene 
expression  changes.  Therefore,  in  subsequent  statistical 
analysis,  we  grouped  datasets  from  the  three  time  points  for 
each  treatment  group  together,  and  treated  them  as  repli¬ 
cates  within  each  treatment  group. 

Hierarchical  clustering  analysis  using  these  genes  showed 
that  study  samples  can  be  grouped  into  three  broad  aggre¬ 
gated  clusters  representative  of  treatment  conditions:  Con¬ 
trol ,  SalGId ,  and  Spor  with  IrradSpor  together,  indicating 
that  gene  expression  changes  caused  by  Spor  and  IrradSpor 
exposure  were  very  similar  (Fig.  2 A).  Gene  expression 
changes,  at  least  until  180  min  post-infection  (the  longest 
time  point  studied),  indicated  that  hepatic  cell  responses  to 
malaria  sporozoite  infection  were  nearly  identical,  whether 
sporozoites  were  irradiated  or  not.  Principal  component 
analysis  of  the  data  (Fig.  2 B)  also  clustered  study  samples 
into  three  well  defined  groups,  with  Spor  and  IrradSpor 
grouping  together.  These  analyses  highlighted  the  close  sim¬ 
ilarity  between  liver  gene  transcriptional  responses  to  Spor 
and  IrradSpor ,  and  their  significant  differences  from  Control 
and  SalGId  groups. 
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FIGURE  3.  Distribution  of  differentially  expressed  genes  across  the  four  study  groups.  A,  number  of  genes  that  are  differentially  expressed  among  the  four 
study  groups.  A  univariate  random  model  F-test  {p  <  0.001 }  was  used  to  determine  the  set  of  genes  differentially  expressed  across  the  four  study  groups.  A  post 
hoc  pairwise  comparison  was  then  used  to  identify  the  subset  of  genes  differentially  regulated  between  different  pairs  of  study  groups.  Each  box  shows  the 
number  of  genes  statistically  different  in  a  pairwise  group  to  group  comparison.  Total  number  of  DEGs  analyzed  (951  genes)  is  shown  in  the  white  boxes. 
Up-regulated  (>2-fold)  genes  in  a  Class  A  versus  Class  B  comparisons  are  shown  in  gray  boxes.  Numbers  in  parentheses  within  each  box  indicates  the  total 
number  of  genes  differentially  regulated  in  each  pairwise  comparison.  8,  scatter  plots  of  pairwise  comparison  among  the  four  study  groups  were  plotted  using 
the  951  DEGs. 


Identification  of  Genes  Associated  with  Sporozoite  Infection — 
Altogether,  a  total  of  951  genes  that  were  differentially 
expressed  among  the  four  treatment  classes  were  identified 
using  a  univariate  random  variance  model  F-test.  A  post  hoc 
pairwise  comparison  was  then  used  to  identify  the  subset  of 
genes  differentially  regulated  between  different  pairs  of  infec¬ 
tion  groups,  and  results  are  summarized  in  Fig.  3 A.  Of  the  426 
genes  altered  in  the  Spor  group  (compared  with  Control ),  107 
genes  were  regulated  with  a  >2-fold  change  in  either  direction, 
with  36  genes  up-regulated  and  71  genes  down-regulated.  In 
the  IrradSpor  group  (compared  with  Control),  426  genes  were 
differentially  expressed,  of  which  24  genes  were  up-regulated 
and  100  genes  down-regulated.  Sporozoite-infected  groups  and 
the  SalGld group  had  more  down-regulated  genes  than  up-reg¬ 
ulated  genes  when  compared  with  the  Control  group.  This 
showed  that  non-sporozoite  containing  salivary  gland 
extract  alone  ( SalGld)  can  induce  gene  expression  changes 
when  added  to  HepG2-A16  cell  cultures  (Scatterplot  of  951 


differentially  expressed  genes  (DEGs)3  for  each  pairwise 
comparison;  Fig.  3 B,  panel  III).  Gene  expression  changes 
induced  by  Spor  and  IrradSpor  treatments  were  almost  iden¬ 
tical  (Fig.  'SB,  panel  V). 

DEGs  from  pairwise  comparisons  identified  regulated  genes 
(relative  to  Control)  that  were  unique  to  each  of  the  three  treat¬ 
ment  groups  {Spor,  IrradSpor,  and  SalGld).  Fig.  4  demonstrates 
that  57  genes  were  uniquely  regulated  in  the  Spor  group  (rela¬ 
tive  to  Control),  compared  with  65  and  232  genes  for  IrradSpor 
and  SalGld,  respectively. 

Pairwise  comparison  between  the  Spor  and  SalGld  groups 
identified  742  selectively  altered  genes  from  wild-type  sporozo¬ 
ite  infection  alone  (supplemental  Table  SI),  244  genes  up-reg- 
ulated  and  227  genes  down-regulated  at  >2-fold  following 
sporozoite  infection  (Fig.  5 A).  DEGs  were  defined  using  Gene 


3  The  abbreviations  used  are:  DEG,  differentially  expressed  gene;  CSP,  circum¬ 
sporozoite  protein;  HSPG,  heparin  sulfate  proteoglycan;  HCV,  hepatitis  C  virus. 
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Ontology  annotations,  to  infer  the  biological  relevance  of  these 
DEGs  to  sporozoite  invasion  of  liver  cells,  and  the  major  Gene 
Ontology  categories  involved  were  those  associated  with 
mitotic  cell  cycle  regulation,  chromatin  packaging  and  remod¬ 
eling,  immunity  and  defense,  protein  proteolysis,  and  intracel¬ 
lular  protein  trafficking.  DEGs  were  also  mapped  to  KEGG  (30) 
and  BioCarta  pathways  (supplementary  Table  S2)  to  under¬ 
stand  their  differential  regulation  within  the  context  of  biolog¬ 
ical  pathways  and  functional  inter-relationships.  The  func¬ 
tional  families  involved  were  those  associated  with  antigen 
processing  and  presentation,  ubiquitin-mediated  proteolysis, 
glycerolipid  metabolism,  extracellular  matrix-receptor  interac¬ 
tion,  cell  communication,  MAPK-signaling  pathway,  fatty  acid 
metabolism,  cell  adhesion  molecules,  xenobiotics  metabolism 
by  cytochrome  P450  and  pathways  associated  with  immune 
system  and  defense,  as  well  as  cell  cycle  and  DNA  replication 
pathways  (Table  1). 

Expressions  of  18  genes  were  highly  differentially  regulated 
(>5-fold  change)  between  Spor  and  SalGld  (Fig.  5 A  and  sup¬ 
plemental  Table  SI).  Three  (NDC80,  SPC25,  NCAPG)  of  the 


Spor  Irrad 


FIGURE  4.  Venn  diagram  showing  co-regulatory  relationships  of  the  DEGs 
among  the  three  study  groups.  DEGs  from  pairwise  comparison  among 
each  of  the  three  study  groups  {Spor,  IrradSpor,  and  SalGld)  relative  to  control 
were  intersected  with  each  other  to  show  the  number  of  common  or  unique 
genes  that  are  regulated. 


four  Spor  up-regulated  genes  encode  nuclear  proteins  involved 
in  the  regulation  of  the  cell  cycle  phase  during  mitosis.  The 
other  up-regulated  gene  (KIF20A)  encodes  for  a  kinesin-like 
motor  protein  that  interacts  with  the  GTP-bound  forms  of 
RAB6A  and  RAB6B  to  regulate  transport  of  Golgi  vesicles  along 
microtubules.  Interestingly,  a  recent  study  on  62  Gabonese 
children  with  varying  severity  of  malaria  (31)  found  that  10 
genes,  including  KIF20A ,  located  in  chromosomal  region 
5q31-33,  were  strongly  correlated  with  malaria  caused  by  par¬ 
asites  expressing  specific  subtypes  of  P.  falciparum  erythrocyte 
surface  protein  1.  Significantly,  KIF20A  expression  was 
increased  and  positively  correlated  with  the  upsA  P.  falciparum 
erythrocyte  surface  protein  1  var ,  gene  subtype. 

Of  the  14  highly  down-regulated  genes,  10  were  for  secreted 
proteins  that  were  either  cytokines/chemokines  (GDF15, 
CTGF,  CCL20,  CXCL2)  or  binding  proteins  (SIOOP,  CP, 
1GFBP1 ,  TTR,  AZGP1 ).  Another  gene,  SERPINE1 ,  encodes  for 
the  serine  protease  inhibitor  (plasminogen  activator  inhibitor 
type  1)  of  tissue  plasminogen  activator  and  urokinase,  and 
functions  principally  in  fibrinolysis,  cell  adhesion,  and  cell 
migration.  In  addition,  plasminogen  activator  inhibitor  type  1 
has  an  important  anti-parasitic  effector  function  and  plasmin¬ 
ogen  activator  inhibitor  type  1 -deficient  mice  showed 
increased  susceptibility  to  P,  chabaudi  blood-stage  malaria 
(32). 

Other  highly  down-regulated  genes  were  CYP1A1,  MAFF, 
TUBB2B,  and  ITIH3.  The  cytochrome  P450  enzymes  are 
monooxygenases  that  catalyze  many  reactions  involved  in  drug 
metabolism  and  the  synthesis  of  cholesterol,  steroids,  and  other 
lipids.  Expression  and  activities  of  some  of  these  liver  enzymes 
are  known  to  be  differentially  regulated  in  malaria.  For  exam¬ 
ple,  cytochrome  P450  family  1,  subfamily  A  enzymes  (CYP1A) 
are  reduced  in  mouse  malaria  models  2-  6  days  following  infec¬ 
tion  corresponding  to  the  period  of  peak  parasitemia  (33, 34).  In 
our  human  liver  infection  model,  down-regulation  of  liver 
CYP1A1  expression  occurred  very  early  (<30  min)  after  sporo¬ 
zoite  infection,  and  is  isoform-specific,  because  expression  of 
CYP1A2  and  CYP2B6  remained  unchanged.  This  discrepancy 


(A)  Spor  Vs  SalGld  (B)  Spor  \/s  Irrad 


FIGURE  5.  Volcano  plot  of  951  DEGs  comparing  the  size  of  the  fold-change  to  statistical  significance  level  in  the  pairwise  comparison  between  (A)  Spor 
to  SalGld,  and  {B)  Spor  to  Irrad .  Vertical  and  horizontal  dotted  lines  indicate  DEGS  with  >2-fold  expression  change  and  p  <  0.001  value  cutoff,  respectively.  A 
random  variance  model  two-sample  t  test  was  used.  Fold-change  was  calculated  for  Spor  relative  to  SalGld  {panel  A)  or  Irrad  {panel  B).  Seven  up-regulated  and 
three  down-regulated  genes  showing  unusually  high  significance  (p  <  10-7)  in  the  Spor  versus  SalGld  comparison  are  annotated  {panel  A).  Twenty  nine  DEGs 
are  found  to  be  down-regulated  in  Spor  versus  Irrad  {panel  B),  and  are  listed  in  Table  2. 
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TABLE  1 

Functional  classification  of  genes  differentially  expressed  genes  between  the  Spor  and  SalGId  groups  using  KEGG  and  Biocarta  pathways 

Gene  expression  changes  were  expressed  as  fold-change  between  Spor  to  SalGId,  and  represent  the  ratio  of  the  geometric  mean  intensities  of  Spor  to  that  of  SalGId .  Datasets 
representing  the  three  time  points  within  each  experimental  group  have  been  grouped  together  and  were  regarded  as  experimental  replicates. 


Pathway 

Pathway  description 

Number  of  genes 

p  value 

1 

hsa04612 

Antigen  processing  and  presentation 

13 

0.000013 

2 

hsa04120 

Ubiquitin  mediated  proteolysis 

5 

0.0000203 

3 

hsa00561 

Glycerolipid  metabolism 

6 

0.0000244 

4 

hsa04110 

Cell  cycle 

26 

0.0001556 

5 

hsa04512 

Extracellular  matrix-receptor  interaction 

7 

0.000186 

6 

hsa01430 

Cell  communication 

6 

0.000171 

7 

hsa04010 

MAPK  signaling  pathway 

13 

0.0000013 

8 

hsa00071 

Fatty  acid  metabolism 

5 

0.000015 

9 

hsa04940 

Type  I  diabetes  mellitus 

5 

0.0004183 

10 

hsa04514 

Cell  adhesion  molecules  (CAMs) 

8 

0.0001145 

11 

hsaOOOlO 

Glycolysis/gluconeogenesis 

8 

0.0000349 

12 

hsa04610 

Complement  and  coagulation  cascades 

11 

0.0001723 

13 

hsa04810 

Regulation  of  actin  cytoskeleton 

7 

0.0000267 

14 

hsa04650 

Natural  killer  cell-mediated  cytotoxicity 

5 

0.0004625 

15 

hsa00240 

Pyrimidine  metabolism 

9 

0.0000175 

16 

hsa04640 

Hematopoietic  cell  lineage 

6 

0.00001 

17 

hsa04060 

Cytokine-cytokine  receptor  interaction 

12 

0.0000817 

18 

hsa00230 

Purine  metabolism 

7 

0.0000133 

19 

hsa04510 

Focal  adhesion 

7 

0.0001511 

20 

hsa04540 

Gap  junction 

5 

0.0002494 

21 

hsa 04350 

TGF-/3  signaling  pathway 

8 

0.0004184 

22 

hsa00564 

Glycerophospholipid  metabolism 

6 

0.0005068 

23 

hsa04310 

Wnt  signaling  pathway 

7 

0.0000095 

24 

h_gl  Pathway 

Cell  cycle:  Gj/S  checkpoint 

6 

0.0002942 

25 

h_mcmPathway 

CDK  regulation  of  DNA  replication 

8 

0.0003691 

26 

h_pparaPathway 

Mechanism  of  gene  regulation  by  peroxisome  proliferators  via  PPARt* 

8 

0.0003646 

may  be  due  to  species  differences  in  malaria  parasite-host 
interactions. 

Tubulin,  the  product  of  the  TUBB2B  gene,  is  an  important 
constituent  of  microtubules  that  act  as  scaffolds  to  determine 
cell  shape  and  as  a  backbone  for  cell  organelles  and  vesicles  to 
move.  ITIH3  is  part  of  the  inter-a-trypsin  inhibitor  family  of 
structurally  related  plasma  serine  protease  inhibitors  involved 
in  extracellular  matrix  stabilization,  and  may  act  as  a  carrier  of 
hyaluronan  in  serum  or  as  a  binding  protein  between  hyaluro- 
nan  and  other  matrix  proteins  to  regulate  the  localization,  syn¬ 
thesis,  and  degradation  of  hyaluronan.  MAFF  is  a  basic  leucine 
zipper  transcription  factor  that  lacks  a  transactivation  domain 
and,  because  it  can  both  homodimerize  with  itself  and  het- 
erodimerize  with  other  leucine-zipper  transcription  factors,  it 
is  a  transcriptional  regulator,  its  role  in  malaria  is  unknown. 

Differential  Effects  of  Wild- type  and  Irradiated  Sporozoites  on 
Hepatocyte  Gene  Expression — A  pairwise  comparison  between 
the  Spor  and  IrradSpor  groups  were  made  using  the  list  of  951 
DEGs  obtained  from  the  initial  analysis  of  variance  analysis 
(Fig.  55),  to  determine  the  subset  of  hepatocyte  genes  whose 
expressions  were  differentially  affected  by  wild-type  sporozo¬ 
ites  compared  with  irradiated  sporozoites.  Twenty-nine  genes 
were  differentially  affected  by  irradiation  of  malarial  sporozo¬ 
ites  (T able  2);  all  were  up-regulated  in  IrradSpor  compared  with 
Spor  (supplemental  Fig.  SI).  Seven  of  those  genes  (. ATAD2 , 
ZNF83,  SMARCA4,  TOPI ,  RBM39 ,  THOC2 ,  and  SART3)  are 
involved  in  nucleic  acid  binding,  suggesting  that  DNA  replica¬ 
tion  and  RNA  transcription/processing  is  the  main  functional 
difference  between  Spor  and  IrradSpor  infections  of  the  liver. 
Other  predominant  biological  processes  represented  by  these 
29  genes  were  macro  molecule  localization  (7  genes),  protein 
transport  and  localization  (10  genes),  and  morphogenesis  of 
anatomical  structures  (6  genes). 
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DISCUSSION 

The  liver  is  the  site  of  development  of  infectious  merozoites 
from  pre-erythrocytic  Plasmodium  sporozoites,  but  little  is 
known  about  the  liver  hepatocyte  responses  in  the  early  stages 
of  infection  with  P.  falciparum  sporozoites.  Our  study  docu¬ 
mented  the  early  gene  expression  profile  changes  with  infection 
of  hepatocytes.  Sporozoite  infection  of  liver  cells  is  rapid  and 
the  in  vitro  inhibition  of  sporozoite  invasion  assay  and  the  inhi¬ 
bition  of  liver  stage  development  assay  using  cultured  hepato¬ 
cytes  normally  require  a  maximum  of  180  min  for  maximal 
sporozoites  binding/invasions  (35,  36).  Therefore,  in  our  stud¬ 
ies,  we  used  a  time  series  experimental  design  to  determine 
hepatic  gene  expression  changes  at  30,  90,  and  180  min  after 
adding  the  P,  falciparum  sporozoites  to  our  HepG2-A16  cul¬ 
tures.  However,  the  different  exposure  times  we  used  had  no 
significant  effects  on  gene  expression  changes.  Hepatic  gene 
transcriptional  responses  were  rapid  and  observed  within  30 
min  following  sporozoite  exposure. 

Many  of  the  differentially  expressed  genes  identified  regulate 
diverse  biological  functions,  such  as  antigen  processing  and 
presentation,  cell  adhesion,  metabolism  pathways,  and  cell 
cycle  regulation.  Although  liver  stage  malaria  is  clinically 
asymptomatic,  a  number  of  sporozoite  ligand-liver  receptor 
interactions  occur  and  enable  sporozoites  to  develop  within  the 
liver  cells  without  initiating  a  local  inflammatory  response  (37). 
Plasmodium  circumsporozoite  protein  (CSP)  is  the  most 
prominent  antigen  covering  the  sporozoite  surface  and  is  crit¬ 
ical  for  sporozoite  function  and  hepatocyte  invasion.  Interac¬ 
tion  of  CSP  with  hepatocytes  inhibits  target  cell  protein  synthe¬ 
sis  (38),  and  alters  the  expression  of  a  large  number  of  genes 
involved  with  various  metabolic  processes,  as  well  as  down¬ 
regulating  genes  controlled  by  NF-k/3  (37).  Our  studies  demon¬ 
strate  that  the  majority  of  hepatic  genes  are  down-regulated 
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TABLE  2 

Liver  genes  whose  expression  is  differentially  altered  upon  exposure  to  radiation-attenuated  sporozoite  ( IrradSpor )  compared  to  wild  type 
sporozoite  (Spor) 

Fold-change  indicated  represents  the  ratio  of  the  geometric  means  oi IrradSpor  to  Spor.  Datasets  representing  the  different  time  points  within  each  experimental  group  have 
been  grouped  together  and  were  regarded  as  experimental  replicates. 


Gene  symbol 

Description 

Parametric  p  value 

Fold-change 

TPR 

Translocated  promoter  region  (to  activated  MET  oncogene) 

1.92e-05 

1.74 

DKFZP434L187 

Hypothetical  LOC26082 

2.75e-05 

1.94 

AP3D1 

Adaptor-related  protein  complex  3,  delta  1  subunit 

8.21e-05 

1.83 

SLC25A36 

Solute  carrier  family  25,  member  36 

8.76e-05 

1.79 

RIF1 

RAP1  interacting  factor  homolog  (yeast) 

9.62e-05 

2.02 

TOPI 

Topoisomerase  (DNA)  1 

0.0001056 

2.46 

SDCCAG10 

Serologically  defined  colon  cancer  antigen  10 

0.000139 

1.8 

SART3 

Squamous  cell  carcinoma  antigen  recognized  by  T  cells  3 

0.0001545 

1.95 

GPC3 

Glypican  3 

0.0001556 

1.8 

Clorfl04 

Chromosome  1  open  reading  frame  104 

0.000156 

1.89 

EHBP1 

EH  domain  binding  protein  1 

0.0001801 

2.27 

RASEF 

RAS  and  EF-hand  domain  containing 

0.0001905 

2.03 

THOC2 

THO  complex  2 

0.0002131 

1.54 

KCNABl 

Potassium  voltage-gated  channel,  shaker-related  subfamily,  (3  member  1 

0.0002981 

1.81 

CDC27 

Cell  division  cycle  27  ho  mo  log  ( Saccharomyces  cerevisiae ) 

0.0003108 

2.17 

NARG1 

NMD  A  receptor  regulated  1 

0.0003196 

1.83 

BAT2D1 

BAT2  domain  containing  1 

0.0004393 

1.77 

MAP4 

Microtubule-associated  protein  4 

0.0004831 

2.46 

PCYOX1 

Prenylcysteine  oxidase  1 

0.0005171 

1.73 

PTPLB 

Protein-tyrosine  phosphatase-like  (prolinc  instead  of  catalytic  arginine),  member  b 

0.0005192 

2.77 

SPG7 

Spastic  paraplegia  7  (pure  and  complicated  autosomal  recessive) 

0.000659 

1.73 

RBM39 

RNA  binding  motif  protein  39 

0.0006645 

2.51 

ATAD2 

ATPase  family,  AAA  domain  containing  2 

0.0006707 

1.54 

HIST1H4C 

Histone  duster  1,  H4c 

0.0006731 

1.58 

SMARCA4 

SWI/SNF-related,  matrix-associated,  actin-dependent  regulator  of  chromatin,  subfamily  a,  member  4 

0.0007271 

1.47 

ZNF83 

Zinc  finger  protein  83 

0.0008767 

2.45 

ANTKRD36B 

Ankyrin  repeat  domain  36B 

0.0008908 

1.88 

ANKRD36 

Ankyrin  repeat  domain  36 

0.0009445 

2.11 

RSRC1 

Arginine/serine-rich  coiled-coil  1 

0.000971 

1.62 

following  sporozoite  invasion.  Similarly,  Singh  etal  (37)  docu¬ 
mented  that  forced  expression  of  CSP  in  HeLa  cells  down-reg¬ 
ulated  10  times  as  many  genes  as  were  up-regulated.  Of  the  742 
differentially  expressed  genes  we  identified,  96  genes  were  sim¬ 
ilar  to  those  reported  by  Singh  et  al  (37)  suggesting  that  the 
altered  expression  of  this  subset  of  genes  can  be  attributable  to 
the  interaction  of  CSP  with  liver  cells  alone,  independent  of  the 
sporozoites  functional  metabolic  processes.  These  genes  were 
predominantly  involved  with  general  responses  to  stimuli  and 
wounding,  protein  localization,  and  intracellular  transport, 
processes  consistent  with  the  role  of  CSP. 

GSTA2,  the  gene  for  glutathione  S-transferase  A2,  was  up- 
regulated  by  exposure  of  HepG2-A16  cells  to  P.  falciparum 
sporozoites  (supplemental  Fig.  S2).  Glutathione  S-transferase 
A2,  a  member  of  the  a  class  glutathione  S-transferases  (GST-a), 
is  the  most  abundantly  expressed  GST-a  in  liver  and  is  released 
rapidly  in  large  quantities  into  the  bloodstream  during  hepato¬ 
cellular  damage  (39).  Besides  metabolizing  bilirubin  and  anti¬ 
cancer  drugs  in  the  liver,  GST-a  exhibit  glutathione  peroxidase 
activity,  which  protects  cells  from  reactive  oxygen  intermedi¬ 
ates  and  peroxidation  products.  Unlike  the  progression  of 
malaria  in  infected  humans  where  circulating  glutathione 
S-transferase  levels  were  decreased  (40),  the  induced  expres¬ 
sion  of  GSTA2  in  sporozoite-infected  HepG2-A16  cells  may  be 
indicative  of  oxidative  stress  contributing  to  the  cells'  inability 
to  support  the  complete  intrahepatic  development  of  P.  falcip¬ 
arum  sporozoites  into  infectious  merozoites.  Because  the  half- 
life  of  GST-a  in  plasma  is  short  (~1  h)  (41),  its  concentration  in 
plasma  will  better  reflect  changes  in  hepatocellular  damage 
than  the  two  widely  used  hepatic  markers,  aspartate  amino¬ 
transferase  (EC  2.6. 1.1)  and  alanine  aminotransferase  (EC 
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2.6. 1.2),  which  have  plasma  half-lives  of  —17  and  47  h,  respec¬ 
tively  (39).  GST-a  may,  therefore,  be  one  suitable  biomarker  for 
early  detection  of  P.  falciparum  infection. 

Radiation  attenuated  whole  sporozoites  have  proven  to  be 
very  effective  in  conferring  sterile  protection  against  malaria 
infection  (42).  Attenuated  sporozoites  are  metabolically  active 
and  their  invasion  into  hepatocytes  and  subsequent  transfor¬ 
mation  into  trophozoites  is  similar  to  that  observed  for  wild- 
type  sporozoites  (43).  However,  unlike  non-irradiated  sporozo¬ 
ites,  trophozoites  from  irradiated  sporozoites  persist  within  the 
cells  and  do  not  transform  into  schizonts  (43).  Although 
HepG2-A16  cells  do  not  support  the  transformation  of  either 
wild- type  or  irradiated  P.  falciparum  sporozoites  into  infec¬ 
tious  exoerythrocytic  forms,  we  did  observe  a  subset  of  genes 
that  were  differentially  expressed  during  the  early  phase  of 
host-pathogen  interactions  between  the  Spor  and  IrradSpor 
groups.  Pairwise  comparison  of  genes  between  the  two  sporo¬ 
zoites  groups  identified  29  differentially  expressed  genes,  all 
up-regulated  in  IrradSpor.  Functionally,  these  genes  were  pre¬ 
dominantly  involved  in  DNA  metabolic  processes,  binding, 
localization,  and  transport  systems  that  are  consistent  with  the 
ability  of  irradiated  sporozoites  to  influence  modifications  of 
hepatocytes  cellular  processes,  which  enable  them  to  reside 
intracellularly  for  longer  periods  of  time. 

Two  major  surface  proteins  of  P.  falciparum  sporozoites, 
CSP  and  sporozoite  surface  protein  2  (SSP2/TRAP),  are  known 
to  have  high  affinity  binding  for  heparin  sulfate  proteoglycans 
(HSPGs)  present  on  hepatocytes  (44)  to  facilitate  the  invasion 
of  sporozoites  into  hepatocytes.  However,  the  class  of  HSPGs  to 
which  they  bind  is  not  known.  Through  binding  with  CSP, 
malarial  sporozoites  exploit  the  sulfation  level  of  host  HSPGs  to 
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navigate  within  the  mammalian  host  (45).  Cells  expressing  low 
sulfated  HSPGs,  such  as  those  in  skin  and  endothelium,  permit 
sporozoite  migration,  whereas  cells  with  highly  sulfated 
HSPGs,  such  as  hepatocytes,  inhibit  migration  and  promote 
invasion  by  sporozoites  (45).  One  important  observation  of  our 
present  study  is  that  the  glypican-3  gene,  GPC3,  is  up-regulated 
in  HepG2-A16  cells  infected  with  irradiated  but  not  wild-type 
sporozoites.  Glypican-3  is  a  heparin  sulfate  proteoglycan  (46) 
secreted  in  the  plasma  of  hepatocellular  carcinoma  patients, 
and  regarded  as  a  diagnostic  serum  marker  for  hepatocellular 
carcinoma  (47-50).  Unlike  wild-type  sporozoites,  irradiated 
sporozoites  are  believed  to  invade  and  persist  within  hepato¬ 
cytes  for  some  time  and  do  not  transform  into  blood  stage  par¬ 
asites.  Up-regulation  of  GPC3  in  irradiated  sporozoite-infected 
cells  could  alter  the  nature  of  the  interaction  between  the 
sporozoites  and  host  cells  to  make  the  hepatocyte  environment 
more  conductive  for  their  long-term  survival  without  further 
transformative  development.  Immunohistopathology  of 
benign  human  liver  tissues  infected  with  hepatitis  C  virus 
(HCV)  showed  strong  staining  for  glypican-3  (51),  and  yeast 
two-hybrid  assays  revealed  that  glypican-3  interacts  with 
CD81,  a  member  of  the  tetraspanin  family  reported  to  be 
involved  in  HCV  infection  and  cell  proliferation  (52).  Indeed, 
hepatocyte  CD81  has  also  been  shown  to  be  required  for  Plas¬ 
modium  sporozoite  invasion  of  hepatocytes  (7,  53),  although 
the  specific  signaling  pathways  for  CD81  interaction  with 
sporozoites  is  unknown.  Recently,  a  prospective  longitudinal 
study  of  a  cohort  of  319  individuals  showed  that  HCV  infection 
leads  to  the  slower  emergence  of  P.  falciparum  parasites  in 
blood  (54).  It  is  possible  that  glypican-3  and  CD81  could  pro¬ 
vide  the  functional  link  between  HCV  and  P.  falciparum  infec¬ 
tions,  and  that  increased  glypican-3  caused  by  HCV  infection 
may  retard  the  transformation  of  P.  falciparum  sporozoites  into 
mature  schizonts.  In  addition,  given  the  current  lack  of  a  suita¬ 
ble  biochemical  marker  for  assessing  “vaccine  take”  for  whole 
sporozoite-based  malaria  vaccines,  glypican-3  is  a  promising 
candidate  for  diagnostic  and  therapeutic  use  in  the  manage¬ 
ment  of  malaria.  The  mechanism  of  protection  conferred  by 
attenuated  whole  sporozoite-based  malaria  vaccines  is  not  yet 
fully  understood.  The  presence  of  liver  stage  parasites  in  the 
hepatocytes  seems  to  be  important  for  maintenance  of  protec¬ 
tion  conferred  by  attenuated  sporozoite  vaccines  (55-57)  via 
CD8+  T  cells.  However,  the  mechanism  by  which  antigen(s)  of 
the  existing  parasites  in  the  hepatocytes  are  processed  and  pre¬ 
sented  to  CD8+  T-cells  is  not  clearly  understood  as  liver  is  not 
particularly  enriched  with  immune  effector  cells  other  than 
Kupffer  cells.  Expression  of  glypican-3  is  known  to  stimulate 
the  recruitment  of  macrophages  into  human  hepatocellular 
carcinoma  tissues  (58,  59).  A  possible  similar  mechanism  of 
recruitment  of  macrophages  by  elevated  glypican-3  during 
sporozoite  infection  of  hepatocytes  may  occur  to  process  and 
present  liver  stage  antigens  to  CD8+  T-cells  to  confer  protec¬ 
tion  from  malaria  with  whole  sporozoite-based  vaccines.  This 
should  be  studied  further  as  a  possible  mechanism  of  protection 
with  attenuated  whole  sporozoite  malaria  vaccines. 

A  subset  of  the  Ras  family  of  genes  (RHOA,  RHOBt  RAB4B , 
RAB22A,  RASD1 )  was  down- regulated  upon  sporozoite  (both 
wild-type  and  irradiated)  infection,  with  only  RASEF  being  up- 
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regulated  (supplemental  Fig.  S3).  The  Ras  family  proteins  are 
small  GTPases  involved  in  signal  transduction  pathways  that 
regulate  cell  growth,  differentiation,  and  survival.  Down-regu¬ 
lation  of  Ras  family  genes  arrests  cell  growth  and  increases  apo¬ 
ptosis.  Up-regulation  of  glypican-3  and  down-regulation  of  Ras 
family  genes  in  the  HepG2-A16  hepatocarcinoma  cells  suggest 
that  P.  falciparum  sporozoite  infection  slows  down  cellular 
growth  and  differentiation  in  hepatocytes.  An  earlier  study 
reported  that  sporozoite  infection  induced  apoptosis  in  hepa¬ 
tocytes  and  that  apoptotic  Plasmodium- infected  hepatocytes 
provide  antigens  to  liver  dendritic  cells  (60).  It  remains  to  be 
determined  if  apoptosis  in  plasmodium  sporozoite-infected 
hepatocytes  is  induced  by  interactions  among  sporozoites,  up- 
regulated  glypican-3,  or  down-regulated  Ras  genes. 

The  HepG2-A16  liver  cell  line  used  in  our  study  is  not  per¬ 
missive  for  functional  development  of  P .  falciparum  sporozo¬ 
ites  into  mature  exoerythrocytic  stage  schizonts,  although 
using  low  growth  density,  this  cell  line  has  been  shown  to  sup¬ 
port  the  exoerythrocytic  stage  development  of  P.  vivax  (54). 
Development  of  exoerythrocytic  stage  forms  in  vitro  has  only 
been  demonstrated  in  the  cytoplasm  of  primary  human  hepa¬ 
tocytes  (61, 62),  HHS-102  hepatoma  cells  (63)  and,  recently,  the 
HC-04  liver  cell  line  (17).  Presently,  there  is  no  human  hepato¬ 
cyte  cell  line  available  in  which  a  high  level  of  sporozoite  infec¬ 
tion  and  formation  of  parasitophorous  vacuole  during  transfor¬ 
mation  of  sporozoites  to  exoerythrocytic  stage  of  P.  falciparum 
have  been  demonstrated.  Because  the  goal  of  our  study  was  to 
identify  early  (<3  h)  changes  in  host  gene  expression  following 
the  initial  infection  by  P .  falciparum  sporozoites  rather  than 
mapping  host  gene  changes  leading  to  the  exoerythrocytic  stage 
development,  the  in  vitro  infection  model  of  HepG2-A16  cells 
was  used.  This  hepatoma  line  has  been  extensively  used  in  P. 
falciparum  sporozoite  invasion  assay  studies  (21-24)  in  which 
P.  falciparum  sporozoites  have  been  shown  to  attach  and 
invade  these  cells  but  do  not  develop  into  mature  exoerythro¬ 
cytic  stage  schizonts.  As  this  cell  line  does  not  support  P.  falcip¬ 
arum  exoerythrocytic  stage  development,  this  study  can  serve 
as  a  starting  point  for  future  comparative  differential  analysis  of 
genome-wide  transcriptomic  changes  between  other  permis¬ 
sive  cells  (such  as  primary  hepatocytes  and  HC-04  cells)  to  help 
identify  genes  and  pathways  that  are  specifically  required  for  P. 
falciparum  transformation  and  development,  including  the 
formation  of  a  parasitophorous  vacuole. 

One  limitation  of  our  study  is  the  small  sample  size,  due  in 
part  to  the  difficulties  in  obtaining  sufficient  amounts  of  fresh 
sporozoites  and  the  limited  amounts  of  RNAs  obtained  from 
infected  cultures.  Although  we  have  not  performed  post-anal¬ 
ysis  experimental  validation  (Le.  using  alternative  assay  tech¬ 
niques  such  as  RT-PCR)  of  our  microarray  results  (which,  when 
performed  on  the  same  samples,  would  only  verify  microarray 
methodologies),  our  results  showed  concordance  with  those 
obtained  in  an  earlier  study  (37).  Of  the  742  genes  whose 
expressions  were  selectively  altered  when  infected  with  wild- 
type  P.  falciparum  sporozoites  ( Spor  versus  SalGld)t  —13% 
were  found  to  be  similar  to  the  subset  of  genes  specifically 
induced  by  malaria  circumsporozoite  protein  (37).  This  agree¬ 
ment  of  genes  from  our  results  to  those  reported  previously  (37) 
not  only  provides  functional  verification  of  our  results  and  val- 
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idates  the  general  performance  of  our  study  but  also  provides 
confidence  in  the  overall  data  obtained  (64). 

In  summary,  using  global  gene  expression  profiling,  we 
have  documented  a  list  of  genes  that  are  differentially  altered 
in  human  liver  cells  following  early  infection  by  P.  falcipa¬ 
rum  sporozoites.  Detailed  analysis  of  the  functional  signifi¬ 
cance  of  these  genes  will  help  provide  further  insights  on  the 
biological  mechanisms  by  which  the  liver  influences  infeo 
tivity  and  transformation  of  sporozoites  into  infectious  exo- 
erythrocytic  forms.  As  the  liver  is  likely  the  only  site  for 
pre-erythrocytic  development  of  malaria  sporozoites,  it 
remains  an  important  gateway  for  invading  sporozoites  that 
can  be  targeted  for  intervention  before  their  development 
into  symptomatic  parasites. 
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